Cosmology with exponential potentials 
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We examine in the context of general relativity the dynamics of a spatially flat Robertson- 
Walker universe filled with a classical minimally coupled scalar field <f) of exponential potential 
V(4>) ~ exp(^fj,(f>) plus pressureless baryonic matter. This system is reduced to a first-order ordinary 
differential equation for ^^(w^) or q(wj,), providing direct evidence on the acceleration/deceleration 
properties of the system. As a consequence, for positive potentials, passage into acceleration not 
at late times is generically a feature of the system for any value of fi, even when the late-times 
attractors are decelerating. Furthermore, the structure formation bound, together with the con- 
straints Q m o ~ 0.25 — 0.3, — 1 < w^o < —0.6 provide, independently of initial conditions and other 
parameters, the necessary condition < fi < 1.6\/8vrGjv, while the less conservative constraint 
— 1 < w^, < —0.93 gives < fi < 0.7\/8tvGn- Special solutions are found to possess intervals of 
acceleration. For the almost cosmological constant case fa —1, the general relation f2^(w^,) is 
obtained. The generic (non-linearized) late-times solution of the system in the plane (u^fi^,) or 
(ti>0,g) is also derived. 
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I. INTRODUCTION 



Scalar-field cosmologies continuously receive through 
time special attention. The interest on scalar fields has 
been stimulated when it has been realized that a scalar 
field might be responsible for inflation of the very early 
universe. If its self-interaction potential energy density 
is sufficiently flat, it can violate the strong energy con- 
dition and drive the universe into an accelerated expan- 
sion. Even if the potential is too steep to drive infla- 
tion, such models still have important cosmological con- 
sequences. Recent observational evidence shows that the 
energy density of the universe today may be dominated 
by a homogeneous component with negative pressure 
(quintessence) Slowly rolling self-interacting scalar 
fields may provide in the late-times evolution of the uni- 
verse a dynamical mechanism to achieve a small effective 
cosmological constant. 

Models with a variety of self-interaction potentials 
have been studied and the cosmological conclusions are 
naturally strongly model dependent. There are attempts 
to base the choice of the potential on some fundamen- 
tal dynamical considerations. A class of potentials that 
is commonly investigated and which arises in a num- 
ber of physical situations has an exponential dependence 
on the scalar field. Such potentials are motivated by 
higher-dimensional supergravity, superstring theory, or 
M-theory that have been compactified to an effective 
four-dimensional theory Also, higher-order gravity 
theories are conformally equivalent to general relativity 
plus a scalar field with potential which asymptotically 
tends to an exponential form 
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Since the first inflationary cosmological models were 
presented, most calculations have been performed us- 
ing one of a variety of approximations. While these 
are usually sufficient and often unavoidable, there is 
considerable interest in scalar field cosmologies that 
can be solved exactly. Exponential potentials have 
been extensively considered in cosmology by various au- 
thors, basically in Robertson- Walker models with no 
other source besides the scalar field Solutions 
with anisotropic/inhomogeneous geometries were ob- 
tained in 0. When a perfect fluid is added, only a few 
solution are known (see references f° r relevant 

works). In the context of ekpyrotic or pre-big bang col- 
lapse, exponential potentials were analyzed in [111 ], while 
non-minimal couplings of scalar fields with exponential 
potentials have been considered in Attention has 

also been directed to the possible effect of such a scalar 
field on the growth of large-scale structure in the uni- 
verse ^jjj . Other recent studies on exponential potentials 
have appeared in [l4| . 

In the present paper, we examine the dynamics of 
a spatially flat Robertson- Walker universe filled with a 
classical minimally coupled scalar field of exponential po- 
tential plus pressureless baryonic matter (luminous and 
nonluminous) non-interacting with the field. Clearly, this 
model cannot be used from the beginning of the uni- 
verse evolution, but only after decoupling of radiation 
and dust. Thus, we do not take into account inflation, 
creation of matter, nucleosynthesis, e.t.c. By performing 
changes of phase-space variables, we reduce the system 
of Einstein equations into a first order nonlinear differ- 
ential equation, which can be expressed as an equation 
for Q, c f > {w^) or q(w^). Solving this equation leads to in- 
tegration of the whole system. Furthermore, studying 
this with respect to its acceleration properties, as well 
as with respect to standard observational constraints, re- 
strictions on the parameter determining the steepness of 
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the potential are obtained. For the almost cosmological 
constant case s» —1, the solution Vl^lxv^) is obtained. 
The late-times behaviour ^(u^) of the system is also 
found. 



II. GENERAL ANALYSIS OF THE MODEL 



We consider the action 



1 



S = —^r / RJ—qd x 
2k 2 ' 



l—g d x 



-g d A x . 



(1) 



with the term L mat corresponding to the perfect fluid, 
and we adopt the Robertson- Walker metric 



ds 2 



-dt l 



a 2 (t)jijdx l dx-i 



(2) 



where jij is a maximally symmetric 3-dimensional met- 
ric, and k = —1,0, 1 parametrizes the spatial curvature. 
The Einstein equations for the system of a perfect fluid 
and a scalar field is equivalent to the following set of 
equations 



p + 3-(p + p) 







V) = 



H 



13(P + P4>) 



(3) 



(4) 



(5) 



where (3 = k 2 /3 = 8itGn/3, a dot means differentiation 
with respect to proper time t, and p^ is the energy density 
of the scalar field defined by 



p<t> = U 2 + v{4>) . 



(6) 



Since the perfect fluid and the scalar field are non- 
interacting, each separately satisfies its own conservation 
equation (J3J, (@J respectively. Equation Q is equiva- 
lent to the classical equation of motion A<j) = V'(<j>) of 
the field (j>, where A is the 4-dimensional Laplacian of 
the metric and a prime denotes differentiation with 
respect to cf>. From Eq. @, it arises that the function 
P(/,(t) is monotonically decreasing (increasing) for an ex- 
panding (contracting) universe. 

From Eqs. J3J , J3J, ©, one can derive the Raychaud- 
huri equation 



H 



3l3 -(p + p + <j> 2 )+ k 
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which is also written in the form 

1 + 3w 

= P\ 6V - -lp A - 

a 



(7) 



(8) 



where w = p/ p characterizes the perfect fluid. 
Defining the quantity 

Q = P 2 + P<t> > ( 9 ) 
we obtain from Eqs. Q, J3J that 

g + G-(g- V + pp) = 0. (10) 

a 

We will specify to a flat universe (k = 0) with dust matter 
(w = 0). Equation (fTTTf) . due to is written for an 
expanding universe as 

^p~^V d = (Vq~P^ + P^) 1/2 (Q-V), 

(11) 

where 4> — ±-\/2 A 
form 



- V. Similarly, Eq. |@} takes the 



p; = T3v^9 ( v ^ 3 P^+P0) 1/2 v / P^ 7 - (12) 
Dividing Eqs. fTTJ) and (fT2")l . we get 

p'<t, _P<t>~V 



Q' 



Q-V 



(13) 



The system of Eqs. ©-© is symmetric under time re- 
versal t — > —t, thus a contracting universe arises from 
an expanding one. Additionally the stable fixed points 
of the system for an expanding universe become unsta- 
ble fixed points in a contracting universe and vice-versa, 
thus we do not discuss the collapsing behaviour further. 
We consider a potential with an exponential-like form 



V{4>) = Voe-i* , 



and define 



R= P -± 



\/Q- P<t> 



(14) 



(15) 



V V 
The system of Eqs. (|12|l . (I13|) is written equivalently as 



180(R- 1) 



(mQ)' = 



R' + pR- 2/i 
2(i?-l) 



(16) 



(17) 



under the conditions ±V (R' - pR) < 0, (R' - pR) 2 > 
18{3R(R — 1). One characteristic of this system is that, 
due to the exponential potential, it does not contain the 
independent variable (f> explicitly. Since the system of 
Eqs. ©-0 is symmetric under the simultaneous change 
p — > — p,, 4> — > —cf), we will consider, without loss of 
generality only p > 0. Combining Eqs. i|ltjfl . I|17ll . we 
can eliminate Q, and obtain 



R' - pR = 



(18) 
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4(iZ - 1)R" - 3i?' 2 - 2//(i? - 3)i?' 

-(R-2) ((18/3 - /^ 2 )i? - 18/3) = . (19) 

Equation (|18f) is equivalent to p^ — —p, which gives from 
Eq. Q a static universe a = constant, and is not of 
further interest here. Making the transformation 



x = -R-i , ^ = ^7 



(20) 



the second order differential equation l|19|l reduces to the 
following first order differential equation in V"(x) : 

= (X - 1)(1 - v^x)^ - 2(X - 2)^ 2 - 3^ , (21) 

where v^ 1 + 1 = 18/3/// 2 . In terms of Xj the previous 
constraints take the form ±Vo(ip~ X~ 1) < 0) 

x-i) 2 >(^- 1 +i)x(x + i) ; 

In terms of the original variables, it is 



X 



2V 



(22) 



and thus Vox > 0- Finding solutions of Eq. (|21|l is the 
starting point for running the equations backwards to 
find solutions of the whole system. Equation (|21|l be- 
longs to the class of Abel equations. There is no known 
systematic way for solving such sort of equations. For the 
same problem with the one examined here, another Abel 
equation, inequivalent to (|2*T|) . was obtained in Ref. @. 
The difference is that there, a transformation of the con- 
figuration variables has been performed. 

Given a solution ip( X ) of Eq. i|21|) . one can find all 
the relevant physical quantities in parametric form. It is 
obvious that x is not in general a global time parameter, 
but it remains a good time parameter, as long as x(t) 
is monotonically increasing or decreasing. So, the scalar 
field 4>(x) is found from Eq. 12211 as 



<t>(x) = <l>i + H 1 ^(x)dx 



(23) 



where <p± is integration constant. Then, the scale factor 
(normalized to unity today) is expressed as 



a\x) 



1 



p Q {v- l + l) X e^W 



(24) 



M^-x-i^-O^+iMx + i) 

with po the present value of /?, and the Hubble parameter 



is 



H 2 ( X ) 



18* ' 



^(x)^-i_ x _!)2 



(25) 



The luminosity distance dj, = (1 + z) J t °dt/a, which is 
directly measured, is given by 

j / \ 1 + 2 

d L (z) = ±— — -x 



V2/ip 



1/3 



rxn 



VVox' 



{V 



x-i 



where the redshift z is defined by 1 + z = 1/a, and xo 
denotes the value of x today. The flatness parameters 
^■m = Pp/H 2 , $l<f, = Pp^/H 2 are found as 



n*(x) = 



(^~ 1 + i)x(x + i) 

(^-!-X-l) 2 



\-n m ( x ). (27) 



Similarly, the energy density and the pressure = 
((f) 2 j '2) — V of the scalar field are given respectively by 
MX) = MX + l)e-^ (x) , P^X) = V (x - l)e"^W. 
Besides these parametric expressions, a physically inter- 
esting point is the determination of the equation of state 
V<t> = Pij){P4>) °f the scalar field, or equivalently of the 
function w < p(p < f l ), where = P(f>/p$ is the state parame- 
ter of the scalar field, given by 



x-i 



X 



1 



W<f, 



(28) 



The inequality Vox > translates to < 1 for Vq > 0, 
and \w$\ > 1 for Vq < 0. Then, it is found that 



P<l>(w<p) 



2V 



-2/ 



■ dw& 



H 1 



W<f, 



(29) 



In this language, Eq. (|21|l becomes a second order dif- 
ferential equation for In p^,(w ( j > ). Concerning the problem 
of acceleration, we can find the decelerating parameter 
q = -ti/aH 2 from Eqs. ®, tf ty . (jj fy as 



q= 2 



3(^ 1 + l)x(l-x) _ 1, 
2(^- 1 -x-l) 2 2 { 



tov+Slf) , (30) 



Since the variables Xi ^ are phase space variables, by 
studying Eq. Q21|). it is straightforward from Eqs. (|30|) to 
get information on the acceleration/deceleration proper- 
ties of the cosmological system. It is seen from the first 
equality of Eq. (|3U|I that acceleration is possible only for 
— 1 < W0 < 0, and thus, it is possible only for Vq > 0. 
Additionally, for Vb > 0, since < < 1, it arises 
from the second equality of Eq. (|30|) that acceleration is 
possible only for —1 < < —1/3. This is an exact 
result, without the usual assumption on the dominance 
of the quintessence component in the energy density, and 
of course, it is only a necessary condition. Finally, the 
relation between x and t is found by 



t(x) 



u ± 



5 M0(x)/2 



dx, (31) 



with t\ integration constant. We observe that the number 
of essential integration constants of the system is three, 
namely po^i, and the one contained in ip(x) horn inte- 
gration of Eq. pifl . This is also predicted from the be- 
ginning, since the system consists of two first-order equa- 
tions ©j and one second-order equation Q, giving 
four integration constants, from where one of them is 
taken out due to the time translation. 
In Ref. @, for V Q > and 



2 H 



H 



(32) 
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the system of Eqs. ©-© was written as a plane- 
autonomous system 



= — 3x 



f 1 ,.2 , 3 n 9 „,2 



y* + -x(l+x'-y') (33) 



xy + -y(l + x 2 - y 2 ) . (34) 



dx 
d\na 

dy 

din a 



It was further found that for any value of fj,, there are 
three kinds of fixed points, characterized by the values of 
(x,y): (i) (±1, 0) , (ii) (0, 0) , (iii) ( M , ^180- /i 2 )/3^3 
for fi 2 < 18(3 and 3y/J!(l, 1)/V2j" for ^ 2 > 9/3. The fixed 
points (i) are unstable nodes or saddle points, the fixed 
point (ii) is a saddle point, while the fixed points (iii) 
are stable points or saddles. More specifically, for any 
value of \i (p? ^ 9(3, 72(3/7) there exists one stable fixed 
point, which is accelerating for /i 2 < 6(3, and decelerating 
for /i 2 > 6(3. For Vq < 0, generic solutions recollapse 
to a singularity after a finite lifetime 10]. The relation 
between the variables x, y and %, ip can be found from 
the previous definitions as 



2 / 2 

X = x jy 



= 1 + x{x- 3^/2(3 / n)y~ 2 . (35) 



III. COSMOLOGICAL IMPLICATIONS 

We are not in position to find the generic solution of 
Eq. H21|) . and then, study from Eq. (|30|1 the accelerating 
properties of the system. However, if / = fi m /f2^ = 
(l/n#) — 1 is the fraction of the flatness parameters, we 
are able to study the accelerating properties of the system 
by deriving from Eqs. i|21[l . I|30|l the following first order 
differential equation for the function f{w t j } ) 



1 



df 



IT 



w^(f + l)dw^ ' 1 ^(iz-i + l)(w + l)(/ + l) 
(l-v- 1 )-(l + v- 1 )w^ > 



, (36) 



where e = sgn(Vo(l — w,p)), and then, 2q(w c / ) ) = 1 + 
3w^/(l + f{w t j } )). For Vo > 0, due to the above discus- 
sion on the fixed points, it is obvious that all the solutions 
of the system start their evolution with deceleration at 
X = +oo (wcj, = 1). Concerning the parameter X; we 
can obtain its time evolution as x — 4'(f J 'P(p ~ 3H(fi)/V = 
±H^/2Vx4 } ~ 1 , thus, for Vo > and <j) < 0, the function 
x{t) (and also the function w,p(t)) is monotonically de- 
creasing. This means that all the solutions emanating 
from the repeller (— 1, 0) have an initial period, charac- 
terized by <j> < 0, where x(t) is monotonically decreasing 
from x = +°° to x = (and also w$ monotonically 
decreases from w,/, = 1 to w^, = — 1). We have plot- 
ted the solutions of Eq. ll.'Uill for these initial time inter- 
vals in Figure 1, shown by the family of "circular" lines, 
where an indicative value of the parameter v = 3/4 has 
been used. The deceleration parameter q(w ( j > ) is shown 



in the same figure by the inclined lines passing through 
the point 1/2 of the vertical axis, providing direct evi- 
dence on the acceleration/deceleration properties of the 
system. Lines of the same indication in Figure 1 cor- 
respond to the same solution, while time increases from 
the right to the left along the curves. Using a first-order 
differential equation, since only one initial value enters, 
there is the merit of representing the whole space of so- 
lutions on a plane diagram. It is seen that there is a 
subclass of solutions which enter from initial deceleration 
to an acceleration era. However, there are still solutions 
with w,)) < —1/3 which remain decelerated. It is easy 
to check from Eq. I|36|l that the structure of Figure 1 
is typical for any value of v. Thus, even for v > 1/2, 
where the late-times attractor is decelerating, there are 
certain classes of solutions which have finite accelerat- 
ing intervals not at late times. Actually, the supernovae 
data favour that the present-day accelerated expansion 
is a recent phenomenon z 9= o ~ 0.5 |L5|. Nucleosynthe- 
sis predictions 0], on the other hand, claim that at 95 
per cent confidence level, it is Q<j,(z — 10 10 ) < 0.045. Be- 
sides this, during galaxy formation epoch (z = 2 — 4), the 
quintessence density parameter is fl^ < 0.05 ^t|- This is 
an additional constraint, applicable in the dust regime, 
which restricts the viable solutions. For instance, for the 
indicative value v = 3/4 of Figure 1, no accelerating solu- 
tion can accommodate this bound, which means that the 
present accelerating era, if possible, has to be located in 
some future time having <j) > 0. For different values of v, 
however, the above bound may match with the present- 
day acceleration, already in the epoch with < 0. 
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FIG. 1: The class of all solutions with Vo > 0, 4> < for an 
indicative value of v = 3/4. Family of "circular" lines repre- 
sent (f2 m /f2^)(w^), while inclined lines passing through the 
point 1/2 of the vertical axis represent q(w^). Lines of the 
same indication correspond to the same solution. The scale 
factor increases from the right to the left along the curves. A 
subclass of solutions show passage from deceleration to accel- 
eration not at late times. Similar portraits occur for all values 
of v. 

In order to pass into the region with <f> > 0, an appro- 
priate time variable, covering a region with both signs 
of 4>, is the parameter 17 m , or equally well f2^, or also / 
(we are going to use /). This is due to that, in general, 
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Q m = 3f3 2 pV(x - 1)/H 3 , and thus, for V > 0, x < 1 
(or — 1 < < 0), the functions fl m (i),f(t) are mono- 
tonically decreasing. Equation (|36[1 can, then, be easily 
inverted, and considered as a first order differential equa- 
tion for w+if), and then, 2q{f) = l+3^(/)/(l+/). The 
solutions w c / > (f) ("lower" curves of the same indication), 
q{f) ("upper" curves) are shown in Figure 2 for the value 
v = 3/4, where time increases from the right to the left 
along the curves. Lines of the same indication in Figure 
2 correspond to the same solution, and also, they are the 
continuations of the corresponding curves of Figure 1. It 




FIG. 2: Class of solutions with Vo > 0, — 1 < w^, < for an 
indicative value of v = 3/4. "Lower" lines of the same in- 
dication represent w</>(/), while "upper" ones represent q(f). 
Lines of the same indication correspond to the same solution, 
and they are the continuations of the corresponding curves of 
Figure 1. The scale factor increases from the right to the left 
along the curves. Similar portraits occur for all values of v. 

is well known that a set of complementary observations 
indicate that the matter energy density f2 m o is no more 
than 0.3 01 j while observations of the higher peaks of 
the cosmic background radiation give tt m o = 0.25 [T^|. 
A conservative estimation for the present quintessence 
equation of state is —1 < w^o < —0.6 It is now 

a simple task, using figures similar to Figures 1 and 2, 
to determine for any particular value of u, if there are 
accelerating solutions satisfying these two constraints on 
J] m o, w^o, together with the previous one from structure 
formation. The result of this analysis is that the above 
bounds are satisfied only for v < 3/4, with agreement 
to previous studies on exponential potentials, based on 
numerical analysis of the full second-order equations of 
the system [2(j ■ This range of v corresponds to a scalar 
field dominated attractor (fl^ = 1), with limiting value 
uiff, = (1 — z/ _1 )/(l + v~ r ). Although there is no ex- 
plicit dependence on the scale factor a or the redshift z 
in Eq. Il3(i[l . and therefore in Figures 1 and 2, however, in 
the previous analysis there is the merit of providing direct 
information on the viability of whole spaces of solutions, 
without particular reference to integration constants or 
other parameters of the system. Additionally, the non- 
dependence of w$ on z is profitable for a quantitative 
treatment, as it is known that changes very rapidly 
between the redshifts of various supernovae data and the 



present. 



A. Approximate solution for w — 1 

The values « — 1 are of particular interest, since 
firstly, it is known pbl l2l| to cover various orders of mag- 
nitude in scale factor during the cosmic evolution, and 
secondly, for a class of potentials larger than those ex- 
amined here, a combined analysis [l9| including all three 
peaks in the Boomerang data of the CMBR spectra, and 
the Type la supernovae, gave as best fit f^o = 0.75, 
— 1 < u)0o < —0.93, thus, behaviour similar to a cosmo- 
logical constant. Considering this restrictive bound for 
and following the previous analysis, we are restricted 
to v < 0.08 as necessary condition for the viability of the 
exponential model. In this case, we are in position to 
find the relation f^iu^). Making the transformation 



t = vx~ 1 > t = ( to )" 3/ V- 1 , 

where r\ = sgn(Vo), Eq. Q21f> takes the form 

dT 2g-r; (l _ qflfo - 
d£ 2^ + 4V? 

For x ~ 0, Eq. H38|) is approximated as 

dT 4 



(37) 



(38) 



(39) 



where 



c 



3ri + 2£ 



3X + 2 



1 5 + We/, 



2a/3^4 2V3(r)x)i 2\/3l-^ 
The solution of Eq. is 



1 + ^0 



ci|T + \/3C| 3 = c 2 



T + 7! c 



(40) 



(41) 



where ci, c 2 are integration constants with |ci| + |c2| > 0. 
In terms of the physical variables w^, fi^, Eq. (|41[) takes 
the form 



ci 



1 + U>0 



±6^2(^ + 1) 



1 + Uty 



= c 2 



7 - W0 



G 




1 + w<$, 



(42) 



which is the first integral between u^,f2^ for ~ — 1. 
This equation represents the "bottom" parts of the lines 
w 4>{f) m Figure 2. Eq. I|42|l is a cubic for 
and can, therefore, be solved for fi^io^). Then, from 
Eq. J37J, (|2*5|) one gets the expression ip(x), and plug- 
ging into Eqs. I|23|l - (|31|l all the relevant quantities can be 
obtained, at least as quadratures. 
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B. Special solutions 



C. Late-times approximation 



For two particular values of v, we can immediately find 
from Eq. I|21[) some special solutions of the system 

i/=l/3 : V = 3/(1-3*) (43) 
v = 9/7 : V = 3/(1-X). (44) 

The cases v = 1/3, 9/7 correspond to ^ 2 = 9/3/2, 81/3/8, 
for which values the general solution of the system we 
discuss has been found in Q (for Vq > 0), and || re- 
spectively. For the solution 1|43|) with Vq > it arises 
that </> > 0. Additionally, we can see that the function 
x{t) (and also w^,(t)) is monotonically increasing with 
< x < 1/3, and that this solution has an eternal accel- 
eration for x > (7— v / 33)/24. Explicitly, for the solution 
(I43|l for both signs of Vq the various quantities are given 
in parametric form from Eqs. H23fl - I|30|) as follows 

=0 1 - Ai - 1 ln(l-3 X ) = 0i-^ 1 ln(2ii^) (45) 

V w$ - 1 / 

a3 = 9p e^ X = Qpoe^ 1 1 ~ «j 
Vb (1-3 X ) 2 4V (1 + 2 W0 ) 2 
9x(l + x) _ 9(1+^) 



(l + 3x) 2 2(2 + i« ) 2 
2V l + 2u; 



e^i (i - W0 )2 
2V Q 



e^i \ \/ 2Vb rv ) 
1 27x(l-x) _ 1 , 27w (l + w ) 



1±4/1 



<7 = o 



2 2(1 + 3 X ) 2 



4(2 + u; ) 2 



(46) 
(47) 
(48) 

(49) 
(50) 



t ± 6e^ / Vox =< ± 3e^ / Vb (!+«;*) 



~9 e M0i + 2 A i2v r (t-t 1 )2 



1 + 2^ 



(51) 



where po,<fii,ti are integration constants. The ± sign of 
Eq. i|4*9"|) is independent of that appeared throughout dis- 
cerning the cases with <j) positive or negative. For Vq > 0, 
since p^ > 0, from Eq. I|48|l it is w<f, < —1/2, and the rel- 
evant sign in Eq. I|49l) is the negative one. In agreement 
with the previous discussion, this solution starts deceler- 
ating with Wtj, ~ — 1, and as increases monotonically 
with time, the universe accelerates resulting to — 
-1/2. Actually, since 8q = 31 - 3(8^ + 3^9 - 80^), 
we have acceleration for any Vl^ > 0.37. We can also see 
that this solution possesses an event horizon [22^. since 
/ dt/a(t) ~ x 1/6 2^i(l/6,5/6, 7/6, 3 X ), and the limit for 
X— »T/3 of this function is finite. 

The solution l|44|) can be easily seen that it does not have 
accelerating intervals, so it is not of special interest here. 



Since equation l|21|l contains phase-space variables, it 
is suitable for approximating the general solution around 
the fixed points of the system. We give for Vq > two 
more equivalent forms of Eq. (|21|l . which will be used 
to obtain the general late-times evolution of the system. 
Defining 

X=-^T > Y = x- 3/ ^-\ (52) 



Eq. I|21(l is written as 
,dY e(2-v)x-l 



Y- 



-Y 



(i-f)x- 



dx 2x 5 / 4 (^X + e) 7/4 4;yx 3/2 (^X + e) 5/2 
where e = 1 (resp. —1) for x > v (resp. x < v )- F° r 

e 



, (53) 



X 



x-i 

Eq. lt?T|) takes the form 
,dY ex-1 



Y 



, Y = X - 3/i ^\ 

(v- l)x-e 



-Y 



(54) 



(55) 



dX 2x 5 / 4 (x + e) 7/4 4z^x 3/2 (x + e) 5/2 
where, as before, e = 1 (resp. —1) for x > 1 (resp. 
X<1)- 

For p 2 < 9/3, Eq. I|53(l is approximated close to the 
fixed point (iii) as 



Y- 



,dY 
da 



2-1/ 



-,Y + u w 0, 



where 



iv + 2 (3i/ + 2)x + 3e 



3(2-1/) L f 3/4 x 1/4 (^X + e) 3/4 
2VT _r ^/3^ + 2 3x + 2 



3(2 - u) \ 1/ 3 / 4 
2VT^17 /3i/ + 2 

3(2 - 1/) V^aTr 



X' 
3x 2 



3/4 



2y : 



a: 3 / 2 y/y 



2y/l^v /3v + 2 



v riv 



3(2-i/)V^/4 (l+u;^)3/4(l-^)i/< 



(56) 

(57) 
(58) 
(59) 
,(60) 



while for /z 2 > 9/3, Eq. Q55J1 is approximated close to the 
fixed point (iii) as 



Y- 



,dY_ 
dv 



v — 1 



Y + v « 0. 



where 



5x + 3e 



l--(5 



l--(5 



1--I5 

j/ 



X 1/4 (x + e) 3/4 
3X + 2 

x 3/4 

3o: 2 + 2t/ 2 \ 

5 + w<j, 



(l+ W0 ) 3 /4(l-^)l/4, 



(61) 

(62) 
(63) 
(64) 
(65) 
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For n 2 < 9(3, the solution of Eq. l|SrJ)l is 



ci \Y + y/l-v, 



,1-1/ 



C2 



y + 



(66) 



Replacing Y from Eq. O, £gs. fgj|, &ecome 
late-times first integrals of the system between the vari- 
ables w</,, (1 + z) 3 /H 2 . 



where c\, C2 are integration constants with |ci| + \q\ > 0. 
Concerning Eq. (JBTJ, its solution for 9/9 < p? < 72(3/7 is 



c 3 \Y + (a - T)v\ a - T = c 4 |y + (<r + r)w| CT+r , (67) 
while, for /i 2 > 72(3/7 the solution is 



In Y 2 + 2avY + v 2 



2a Y + av 

z arctan ■ 



VT^o 2 



(68) 

In the above two equations, 03,04, c are integration con- 
stants with I C3 1 + I C4 1 > 0, and a = l/2yl — v~ x , 
t = V(4-3i/)/4(i/-l). 

We can write Y from Eq. I|35[) in terms of x, y as 



y = 



(69) 



and thus, replacing y it, u in terms of x, y from Eq. (16911 . 
Q, respectively, £ 9 s. (EJ), become 

late-times first integrals of the system in the plane 
(4>/H, y/V /H), approaching respectively, for /x 2 < 9(3 a 
scalar field dominated attractor (node), for 9(3 < /j 2 < 
72(3/7 a node "scaling" attractor, and for /i 2 > 72(3/7 a 
spiral "scaling" attractor. 

Furthermore, we can write Y from Eq. 1271) . (I30|l in terms 
of Wcf,, Q,p or u>0, q as 



y = 




(l-^) 1 /MVi + ^ 



Thus, replacing Y, u, v in terms of w^, Q<f, from Eqs. Q7U[I. 
(|5U|l . J55J) respectively, i?gs. fiTS)) . f37| ), f??5)) become late- 
times first integrals of the system between the variables 

Finally, can be written from its definition as Q<f, — 
1 — (3po(l + z) 3 /H 2 , and thus, Y is written as 



y= 



V2 



(1-^) 1/4 IV 1+ ^ V l-/3po(l + ^) 3 ff- 



IV. CONCLUSIONS 



In conclusion, we have studied how long a scalar field 
of exponential potential can play the role of quintessence 
in a realistic universe filled with pressureless baryonic 
matter. By deriving and studying directly the first or- 
der differential equation of the deceleration parameter q 
with respect to the state parameter of the scalar field 
Wff,, we found that there exist whole classes of solutions 
which enter into accelerating eras not at late times, inde- 
pendently of the parameter fj, determining the steepness 
of the potential. Moreover, the condition < —1/3 
does not necessarily imply acceleration. Observational 
constraints of fi m o,iU0o provide immediate bounds on /1 
for the viable solutions, irrespectively of initial condi- 
tions or other parameters, in accordance - as it has been 
claimed |23] - with all the current observations. Using 
the same first order equation, we have obtained the so- 
lution Q^Wfj,) for the almost cosmological constant case 
W(j, w — 1, and the generic late-times solution of the sys- 
tem in the plane of the variables (4>/H, y/V / H), (w^, fi^), 
(w^q), and [wj,, (1 + zf/H 2 ). 



(72) 
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